Transport through molecular junctions with a nonequilibrium phonon population 
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The calculation of the nonlinear conductance of a single-molecule junction is revisited. The self 
energy on the junction resulting from the electron-phonon interaction has at low temperatures log- 
arithmic singularities (in the real part) and discontinuities (in the imaginary one) at the frequencies 
corresponding to the opening of the inelastic channels. These singularities generate discontinuities 
and logarithmic divergences (as a function of the bias voltage) in the low-temperature differential 
conductance around the inelastic thresholds. The self energy also depends on the population of the 
vibrational modes. The case of a vibrating free junction (not coupled to a thermal bath), where 
the phonon population is determined by the bias voltage is examined. We compare the resulting 
zero-temperature differential conductance with the one obtained for equilibrated phonons, and find 
that the difference is larger the larger is the bare transmission of the junction and the product of 
the electron dwell time on the junction with the phonon frequency. 
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Electrons passing through a small molecule can change 
its quantum state and this usually requires a finite en- 
ergy transfer from the transport electron, yielding inter- 
esting structures in the I-V characteristics. This rich 
pattern depends on featuress such as the equilibration 
time of the vibrations compared to the typical time be- 
tween consecutive electrons passing through the junction, 
or whether the electrons can pump more and more exci- 
tations into the vibrational states. Single-molecule junc- 
tions based on direct bonding of a small molecule between 
two metallic electrodes seem by now rather established 
experimentally. There are also quite a number of the- 
oretical studies, focusing on various regimes of the rele- 
vant parameters ; 12-27 In particular, the modification of 
the differential conductance at the opening of the inelas- 
tic channel has been an issue of considerable interest (an 
extended discussion may be found in Ref. [HI). 

An electron crossing the molecular bridge can do so 
with or without changing the vibrational excitation state 
of the molecule. At low temperatures, the first inelas- 
tic channel comes in when the bias voltage, V, exceeds 
huj Q /e, where wo is the normal frequency of the junction. 
This however does not necessarily imply an increase of 
the total conductance, since the elastic conduction chan- 
nel is modified as well. Technically, the low-temperature 
conductance associated with the opening of the inelastic 
channel stems from two sources. The first is the imag- 
inary part of the junction self energy, induced by the 
interaction with the oscillator. This function develops 
discontinuities, so that there appear additional contribu- 
tions to the conductance only for eV > Hloq. Discon- 
tinuities in the imaginary part of the self energy imply 
logarithmic divergences in its real part, via the Kramers- 
Kronig relations ; 25 ' 28 As a result, the low-temperature 



differential conductance develops logarithmic singulari- 
ties (as a function of the bias voltage) around the inelas- 
tic thresholds ; 14 ' 17 ' 25 The second contribution is due to 
a nonequilibrium population of the vibrations. At zero 
temperature, phonons can be excited only when the en- 
ergy of the electrons (i.e., eV) exceeds ujq, leading to an 
additional modification of the conductance at these volt- 
ages. 

In a recent publication-^ we have presented a detailed 
calculation of the differential conductance and analyzed 
its behavior at low temperatures. Our calculation implic- 
itly assumed that the molecular junction is in a good con- 
tact with a thermal bath, such that the vibrations follow 
the Bosc-Einstcin distribution, with the same tempera- 
ture as the two leads. In particular we have considered 
the conditions for the conductance to increase or decrease 
at the channel opening. In this short communication we 
repeat that analysis for the case of a free junction (i.e., 
the oscillator is not coupled to a heat bath), for which at 
low temperatures the phonon population is determined 
by the bias voltage. This case is analogous to the mag- 
netization of a Kondo ion out of equilibrium^ 

In this note we use the same Hamiltonian and the 
same notations as in our previous paper ^ To lowest or- 
der in the electron-phonon interaction on the dot, 7, the 
scattering-down rate (de-exciting the vibration) can be 
derived from the golden-rule, 
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and the scattering- up rate, w up , is obtained from Eq. ([!]) 
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by swapping u>q and — ujq- Here, f a (oj) = (exp[/3(u; — 
+ 1) i s the Fermi distribution in the left (a = L) 
or the right (a = R) lead, in which the chemical potential 
is n a , and N(uj) is the bare density of states on the dot, 
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It is assumed that the dot, modeled by a single level of 
energy e, is coupled symmetrically to the two leads, with 
To being the bare resonance width and eV — — [ir, 
with fiL > /J.R- All electronic frequencies (using h = 
1) are measured from e — /i, where // = (/!£, + [ar)/2. 
The kinetic equation for the vibration population, -/V no , 
is then (see also Ref. QT 



dN ho (t) 
dt 



-N ho (t)w down + [1 + N ho (t)]w 



(3) 



The scattering-up rate is due to all processes by which an 
electron can excite an oscillator mode, losing the energy 
luq in the process, and moving over from the left (right) 
lead back into the left (right) lead, or from the left (right) 
lead into the right (left) one. Hence, at zero tempera- 
ture, w up ^ only when the bias voltage exceeds the 
frequency ujq, resulting from the L — » R process. On the 
other hand, Wdown results from all four processes in which 
the electron gains the energy luq. However, at zero tem- 
perature and when eV > u>o the R — ► L process ceases to 
contribute. It follows from Eq. ([3]) that the stationary 
vibration population, 
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is independent of the electron-phonon coupling^- and 
vanishes at zero temperature as long as eV < loq- 
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FIG. 1: The population of the oscillator modes, A^o, at zero 
temperature, as a function of eV/coo, for various values of the 
bare transmission: T—0.2 (smallest dash size), 4/17, 0.3, 4/9, 
and 0.9 (largest dash size). Here c^o/Fo = 1/3. 

Figures HJ [21 and |3] depict the nonequilibrium phonon 
population at zero temperature (see also Ref. [l3). The 
curves are for various values of the bare transmission, 
T = Tq/ (fi 2 + Tq), of the junction, The population usu- 
ally increases with T and with V, and its magnitude in- 
creases with the ratio ujq/Tq^ This can be understood 



qualitatively: since eV — loq is the driving force pushing 
the population out of equilibrium, A^ increases with the 
voltage. The increase with T is due to the fact that the 
dwell time on the dot, r^, increases with T from being 
very short off resonance to becomings 1 /Tq around the 
resonance. As emphasized in Ref. [2|| to effectively ex- 
cite the phonon, t ( i should be longer than the response 
time of the oscillator (about uJq 1 ), i.e., Tq < ojq. Another 
interesting issue is that the time interval between succes- 
sive electrons passing the junction, t c ~ e/I ~ l/(eVT), 
decreases as T or eV are enhanced. This too will cause 
eV to affect more significantly the population at higher 
values of the bare transmission. However, when t c < t^, 
Pauli constraints between consecutive electrons on the 
dot should come into play. 
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FIG. 2: The same as in Fig. [U for oj = T . 
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FIG. 3: The same as in Fig. [T] for aj /r = 3. 

When r 3> wq, one may estimate the vibration pop- 
ulation at finite temperatures by assuming that the bare 
density of states, [see Eq. ((5J)], is a constant. Then the 
stationary population becomes 

4" 
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independent of the bare resonance width Tq and of the 
band width of the leads. This resembles the magnetiza- 
tion of a Kondo ion out of equilibrium^ with the oscilla- 
tor frequency wo playing the role of the applied magnetic 
field. (That magnetization is the inverse of 1 + 2A?ho-) 
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The current flowing through the junction can be pre- 
sented in the form [see, e.g., Ref. [281 ] 



I = eT r 



did 



7T \ 



(6) 



where ImG^/ir is the density of states on the dot, fully 
dressed by the interaction with the vibrations. The (ad- 
vanced) Green function on the junction is 
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where Aeo is the polaron shift, 



(7) 



Ae = - 
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^(/iHl/sM), (8) 



and Sho is the self-energy due to the electron-phonon 
processes. Up to second-order in the electron-phonon 
coupling it reada2£ 



2/ \ cj — cj n — cj' — i0 + 
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Employing Eqs. ([5]) and © in Eq. ©, we have derived 
the differential conductance of the bridge, G, at zero tem- 
perature, and to second order in the coupling with the 
phonons. 
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FIG. 4: The zero-temperature differential conductance, as 
a function of eV/uo, for two values of the bare transmis- 
sion, marked on the figure. The full lines show the conduc- 
tance with the nonequilibrium phonon population, the dashed 
curves are with the equilibrium one. Here ojo/Fq = 1/3. 

FigureslU [5J and[6]show several examples of the depen- 
dence of the differential conductance on the bias voltage 
and on the other parameters of the junction. The first fig- 
ure corresponds to the case where the junction is tightly 
bound to the leads, and hence the dwell time of the elec- 
trons is rather short. Then, at relatively low values of 
the bare transmission (T = 0.3 in our example) there is 
no discernible modification in the conductance which is 



about the same either for nonequilibrium phonon pop- 
ulation or for the equilibrium one. However, for higher 
values of the bare transmission (T = 0.81) the step-down 
feature of the conductance at threshold for inelastic tun- 
neling is enhanced for nonequilibrium phonons as com- 
pared to the equilibrated ones. 
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FIG. 5: The same as in Fig. 0] for lu =T . 

From the few examples portrayed in Figs. |4]|6]we see 
that when the bare transmission of the junction is high, 
then the conductance in the presence of nonequilibrium 
phonons is lower than the one pertaining to the case of 
equilibrated vibrations. For low bare transmissions the 
difference is rather small. One notes that the logarith- 
mic singularity associated with the real part of the self 
energy at the channel opening, which has been previ- 
ously discusse d 14 ^ 25 ' 28 for an equilibrium population of 
the phonons, is still manifested also when these excita- 
tions are out of equilibrium, for high enough values of 
the bare transmission of the junction. 
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FIG. 6: The same as in Fig. H for m = 3r . 

The fact that the changes in the differential conduc- 
tance which are associated with the type of the vibra- 
tion population, ./Vho, are more pronounced for the larger 
wo /To ratio is connected with the actual value of -/Vho, 
see Figs. QJ3] The electrons pump more and more ex- 
citations into the higher vibrational states as the bias 
voltage increases and this pumping is more effective as 
the dwell time exceeds considerably the response time of 
the oscillator. 
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